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Abstract In this paper, the axion contribution to the electromagnetic wave propagation is studied. 
First we show how the axion electrodynamics model can be embedded into a premetric formalism of 
Maxwell electrodynamics. In this formalism, the axion field is not an arbitrary added Chern-Simon term 
O |' of the Lagrangian, but emerges in a natural way as an irreducible part of a general constitutive tensor. 

We show that in order to represent the axion contribution to the wave propagation it is necessary 
to go beyond the geometric approximation, which is usually used in the premetric formalism. We 
derive a covariant dispersion relation for the axion modified electrodynamics. The wave propagation 
in this model is studied for an axion field with timelike, spacelike and null derivative covectors. The 
birefringence effect emerges in all these classes as a signal of Lorentz violation. This effect is however 
completely different from the ordinary birefringence appearing in classical optics and in premetric 
electrodynamics. The axion field does not simple double the ordinary light cone structure. In fact, 
it modifies the global topological structure of light cones surfaces. In CFJ-electrodynamics, such a 
modification results in violation of causality. In addition, the optical metrics in axion electrodynamics 
are not pseudo-Riemannian. In fact, for all types of the axion field, they are even non-Finslerian. 

C ■ . Keywords Axion electrodynamics • wave propagation • birefringence 
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1 Introduction 



Axion electrodynamics, i.e., the standard electrodynamics modified by an additional axion field [T], is a 
subject of a growing theoretical and experimental interest. In particular, such a modification provides 
a theoretical framework for a possible violation of parity and Lorentz invariance — the Carroll-Ficld- 
Jackiw (CFJ) model [2], [3]. For recent developments of this model, see [4] and the reference given 
therein. Also the non-abelian extensions of the axion modified electrodynamics for the Standard Model 
[5] and gravity [6] were worked out. 

Although Lorentz invariance is a basic assumption of the standard classical and quantum field 
theory, in quantum gravity and string theory this invariance is probably violated. One believes that the 
low-energy manifestation of Lorentz violation can be observed in experiments with the electromagnetic 
waves. It justifies the importance of examination of the theoretical aspects of axion contributions to 
the light propagation effects. Although this problem was investigated intensively, we apply here an 
alternative approach based on a metric-free (premetric) formulation of electrodynamics, see [7] and 
the reference given therein. A characteristic feature of such a construction is that the axion field is not 
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involved by hand (merely as an additional term in the Lagrangian). Alternatively, in the premetric 
formalism, axion emerges as a necessary and natural partner of the standard photon variable. 

In fact, there is a certain contradiction between the premetric electrodynamics and the CFJ-model. 
On one hand, in the premetric construction, one usually states that the axion component does not 
alternate the wave propagation at all [7]. On the other hand, the Lorentz violation in the CFJ-model is 
explicitly manifested in a modification of the standard light cone. In fact, this contradiction comes from 
the specific geometric optics approximation which is usually applied in the premetric electrodynamics. 
A constant axion field indeed does not contribute to the wave propagation. To describe the wave 
propagation in the CFJ-model one has to go beyond the geometric optics approximation and take into 
account the first order derivatives of the axion field [8] , [9] . With this modification, two constructions 
yield the equivalent results. 

The original CFJ-model is based on the standard special relativistic electrodynamics Lagrangian 
modified by an additional Chern-Simon term [2] 

L = -\f^ - ^p^AjdkAi + Af , (1) 

where pi is a constant time-like vector. In fact, a special parametrization pi — (/x, 0,0,0) is usually 
used. Note that, in such a construction, the Lorentz violation is involved by hand. The spatial 50(3)- 
invariance and gauge invariance are preserved however. It was already pointed out in [5] that the model 
can be equivalently reformulated in an explicitly gauge invariant form 

L = ^ /•;.,/• '•• + ^ M F tJ F kl + A a l , (2) 
where ip is a pseudoscalar axion field. It is related to the vector pi as 

= Pi, or ip = lit. (3) 

In the current paper, we use a Lagrangian similar to ([2|). We will not restrict, however, to an 
axion field with constant first order derivatives and will not require it to be timelike. Moreover, we 
will consider a model on a curved pseudo-Riemannian spacetime. Our approximation will be based, 
however, on an assumption that the gravitational field changes much more slowly than the axion field. 
This restriction considerably simplifies the calculations and, hopefully, does not change the results, at 
least, qualitatively. 

The current paper is devoted to the 60th birthday of Professor Bahram Mashhoon. His permanent 
interest and considerable contribution to the study of wave propagation effects are well known, see, 
for instance, [12], [13] . The methods developed by Mashhoon will be useful also in the axion modified 
electrodynamics. 



2 The axion modified electrodynamics model 

Let us briefly describe how the axion modified electrodynamics is embedded in the framework of the 
premetric approach. We start with two independent antisymmetric tensor fields, the field strength 
tensor F^ and the excitation field The latter field is a pseudotensor density of weight (+1). The 
flux conservation law (the first Maxwell equation) is postulated, 

e ijkl F ij>k = 0. (4) 

Here the Roman indices range from to 3, e 1 ^ 1 and eyjw are the Levi-Civita permutation tensors 
normalized with e 0123 = —60123 = 1, the commas denote the ordinary partial derivatives. Due to ([4]), 
the field strength tensor is expressed in term of the standard vector potentials Aj, 

Fii = \ (Aij - . (5) 
A local linear homogeneous constitutive relation between the fields Fy and H 1J , 

U H = I X ««fjy , (6) 
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is assumed. By definition, the constitutive pseudotensor \^ kl must respect the symmetries of the fields 
Fij and H ij , 

yijkl _ y\ij]kl _ ^ij[kl] ^ /ys 

Hence it has, in general, 36 independent components. The irreducible decomposition of this tensor 
under the group of linear transformations involves three independent pieces. One of these three pieces 
is the axion field, which is a subject of our interest. 

The high number of components of allows to describe a wide range of physical effects. The 
axion field, however, adds only one addition component to the standard electrodynamics. So, in order 
to embed the axion electrodynamics into the prcmctric approach, one has to restrict the number of the 
independent components. The first restriction comes from the Lagrangian formulation of the model. 
We assume the action integral to be of the standard form 

A = f (/,/«'•• • A,J'),l'.r. (8) 

J M 

When (|6|) is substituted, the action takes the form 

A = J (he^FijFu + AiJ^j d^x . (9) 

Note, however, that, in contrast to the ordinary textbooks formulation, the electromagnetic current 
J 1 and the excitation field TC 1 ^ are pseudotensorial densities of weight +1 (see Appendix for definition 
and details). Thus the Lagrangian (|8|) is general covariant and admits arbitrary transformations of 
coordinates. The constitutive pseudotensor involved in jSJ respects an additional symmetry 

x ijkl = x klt3 ■ (10) 

This condition removes 15 independent components of x ij which compose the so-called skewon part 
of the constitutive tensor. The remaining quantity of 21 independent components is decomposed irre- 
ducibly to a principle part of 20 components plus one component that represents the axion field. In 
contrast to the Lagrangian ([2]), the premetric formulation ([8]) does not involve an addition axion term 
in the Lagrangian. In fact, such a term is hidden in the constitutive tensor. 

The variation of (JSJ) with respect to the vector potentials Ai yields the second Maxwell equation 

W i J = J i , or l(x ijkl Fu) ij = J i . (11) 

Observe that this general construction, is explicitly gauge invariant. As usually, the charge conser- 
vation law diJ % = is a consequence of the field equation (fTTjl . 

The standard electrodynamics is reinstated in this formalism by a special Maxwell constitutive 
tensor 

(H«) X ii« = ^ik^l _ gtigjk) x — (l2) 

Here g % $ is a metric tensor with the signature {+, — , — , — } and with the determinant g. To involve the 
axion field contribution to the standard electromagnetic Lagrangian, it is enough now to consider a 
slightly modified constitutive tensor of the following form 

^ijkl _ (Max)^rjTd _|_ ^p^ijkl _ 

Here the axion ip is a pseudoscalar field. It is invariant under transformations of coordinates with 
positive determinant and changes its sign under transformations with negative determinant. 

The constitutive tensor () 13[) is not merely a modification of the standard expression (I12p . In fact, 
it is an irreducible decomposition of a general constitutive tensor \ n M ■ The Lagrangian formulation 
removes the skewon part of it. An additional requirement of closeness, see [7], restricts the principle 
part of 20 independent components to a pure metrical expression. So the axion part appears in this 
formalism as a natural and necessary ingredient of a general construction. To simplify our consideration 
we will treat the axion field only phenomenologically. An additional dynamical axion Lagrangian can 
be readily added to ©■ 
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3 Wave solution 

To study the wave propagation in the axion modified model, we assume the electromagnetic current 
to be equal to zero. Substitute (J5]) into the second field equation (fTTjl to rewrite it as 

( X iikl A k ,i) tj = 0, (14) 

or, equivalently, as 

X tfW ^*,ii+X <iW ^,i = 0. (15) 
In this paper, we will apply the following approximation 

X m ,j=t^ m - (16) 

In other words, we restrict to the spacetime regions where the gravitational field varies slowly as 
compared to the change of the pseudoscalar field. In particular, our analysis will be exact for an axion 
field considered on a flat Minkowski space. 
Substituting flT3J) and JTSJ) into (JTSJ) we get 

(9 ik 9 jl - 9 ll 9 jk ) A k ,ijV^9 + ti<? jklA k,l = . (17) 
We are looking for a plane monochromatic wave solution of the equation (|17|) . Write it as 

A k =Rc(a k e l ^ xm ) . (18) 



Here the amplitude a k and the wave covector q m do not depend on a point. Both quantities can be 
complex, the physical solution A k is equal to the real part of the corresponding complex expression. 
Since (fT"7|) is a linear field equation with real coefficients, it is possible to deal, as usual, with the 
complex valued expression A k = ake lqmX itself. Substituting this ansatz in (jTTJ) we have 

(»V - 9 U 9 jk ) V=9Wia k - ii>,W ijkl a k = . (19) 

This tensorial equation is represented by a linear system of four equation for four components of the 
covector a k - Write it briefly as 

M ij aj = Q, (20) 

where 

= (9 ij q 2 ~ </V) ^9 + i^kil^ ■ (21) 

Observe two evident relations that hold due to the definition of the matrix M lk 

M ij qi = 0, M ij qj = 0. (22) 

These relations can be given a pure physical sense: The former one represents the charge conservation 
law, while the latter relation represents the gauge invariance of the field equation. 

The linear system (|20[) has a non-zero solution if and only if its determinant equal to zero. For 
the system (|2U|) . this condition holds identically, which can be seen even without explicit calculation 
of the determinant. Indeed, the identities (I22|) express linear relation between the rows (and between 
the columns) of the matrix M lk . So this matrix is singular. Moreover, ([22l ) also means that the linear 
system (|20l) has a non-zero solution of the form 

aj = Cqj (23) 

with an arbitrary constant C. This solution is evidently unphysical since it corresponds to the gauge 
invariance of the field equations. To describe an observable physically meaningful situation we must 
have an additional linear independent solution of (|2"U1) . 
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4 Dispersion relation 

The linear system (|19j) has two linear independent solutions (one for gauge and one for physics) if and 
only if its matrix M lJ is of rank 2 or less. An algebraic expression of this requirement is 

A, = , (24) 

where Ay is the adjoint matrix. This matrix is obtained by removing the i-th row and the i-th column 
from the original matrix. The determinants of the retaining 3x3 matrices are calculated and assembled 
in a new matrix Ay-. The entries of the adjoint matrix are expressed via the entries of the matrix M %3 

as 

Ay = ^e Ull2l3 e nij2j3 M^M^M^ . (25) 

Note that for a covariant tensor of a rank (2, 0), the adjoint matrix is a contravariant tensor of a rank 
(0, 2). Since the adjoint matrix has, in general, 16 independent components it seems that we have to 
require 16 independent conditions. The following algebraic fact |15) shows that the situation is rather 
simpler. 

Proposition: If a square n x n matrix M lJ satisfies the relations 

M ij qi =0, M ij qj =0 (26) 
for some nonzero vector qi, its adjoint matrix Aij is represented by 

Aj = KlH<lj ■ (27) 



For a formal proof of this fact, see [19j. Consequently, instead of 16 conditions (|24j) . we have to require 
only one condition 

X(q) = 0. (28) 

Recall that this condition is necessary for existence a physically meaningful solution of the generalized 
wave equation, so it is a generalized dispersion relation. 

We calculate now the adjoint matrix for the axion modified electrodynamics model. Substituting 
(|2"T) into ([25]) we get 

Aj = 7y e iiii2i3 e jjlj2j3 (-M ljl + ilp.killl^ ^ 1 

x (M l2n + i^M^ 2nk2h ) Of + it/j M q h e l3: > 3k3h ) , (29) 

where a notation 

M ij = (g^q 2 - q l ql) V=5 (30) 

for the pure Maxwell part of the matrix M y is involved. 

Calculate now in turn the entries of (|2"9"|) as the powers of the imaginary unit. The calculations are 
considerable simplified when we take into account that due to the Proposition above the result has to 
be a symmetric matrix. It means that all antisymmetric contributions to Ay are canceled. The free in 
i term, i.e. the pure Maxwell term, takes the form 

j^i^ihhhM^M^'M^' = -V~gq%qj ■ (31) 
The i-term takes the form 

i(l/2)e« li2i3 €^ 1 ^ 3 M il ^M^^, fe 3* 3 ^' 3fe3 ' 3 - (32) 

Since the Maxwell matrix M lJ is symmetric, this expression is antisymmetric in the indices i,j and 
does not give a contribution to the adjoint matrix. In fact, explicit calculations show that the expression 
(|32|) vanishes identically. The z 2 -term is given by 

i 2 (l/2)e iili2i3 e J - ilj2j3 M ilJ ^ )fe2 « 2 e i2i3fe3!2 ^,fe3%e i3i3fc3i3 
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(</></>< m )^- {ip, m q m f (33) 
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The z 3 -term takes the form 

« 3 (l/2K ni2l3ejJ1 , 2J 3^^^ (34) 

This expression is evidently antisymmetric in the indices i,j so it does not give a contribution to the 
symmetric matrix Aij. In fact, it is equal to zero. 

Consequently we have the adjoint matrix in the following form 

q A + 0>, m ^ m ) q 2 ~ (i>, m q m ) 2 ] QiQj ■ (35) 

Thus the dispersion relation for the electromagnetic waves in the axion electrodynamics is expressed 
as 

9 4 + (M' m ) q 2 ~ (ip,mq m ) 2 = . (36) 

Observe some characteristic features of this equation. 

1. In Minkowski space, for an axion field with a constant covector of derivatives, it coincides with 
the dispersion relation expression given in [2]. 

2. The axion dispersion relation (I36p is essentially different from the general covariant dispersion 
relation appearing in the premetric electrodynamics. The premetric dispersion relation is a quartic 
homogeneous polynomial in the wave covector variable q. Its general form is 

g l3kl q t q 3 q k qi=0. (37) 

Certainly the homogeneity is originated in the geometric approximation used for its derivation. The 
quartic polynomial of the axion electrodynamics ([36| is not homogeneous, so it provides some additional 
types of light cones structure. In particular, the birefringence effect here comes from the derivatives 
of the media parameters. On the contrary, in the premetric electrodynamics as well as in the classical 
crystal optics, the birefringence effect comes from the tensorial nature of the parameters of the media. 

3. Another interesting feature of the dispersion relation p6|) is that it is, in fact, a general covariant 
expression. Indeed it is invariant under general pointwise transformations of coordinates (with positive 
or negative determinant). This is in spite of the fact that we used for its derivation the ansatz (1181) 
which is only special relativistic. 

4. The relation (|36[) is invariant under the transformation q l — > —q % . Thus the light cones structure 
has to be PT invariant. 

5. The relation (136|) is invariant under the transformation of the field ip* — > — ip i . Thus the light 
cones have to be similarly oriented relative to the vector ip l . 

6. A quartic expression dispersion relation (|36[) is sometimes factored to a product of two second 
order polynomials. It is easy to see that the relation (|36[) cannot be factored in covariant way. Such 
factorization is possible, however, in special coordinates. Since (|36|) is not homogeneous, at least one 
of these quadratic factors is not homogeneous. 

7. For ordinary electrodynamics in vacuum, the factors are homogeneous and coincide. This unique 
homogeneous quadratic factor corresponds to unique Minkowski metric. For general anisotropic media, 
in (geometric approximation) two factors are homogeneous but different one from another. This case 
is formulated in term of two optical metrics, both pseudo-Riemannian. In axion electrodynamics, at 
least one factor is necessary inhomogeneous. This factor cannot be reformulated in term of on ordinary 
pseudo-Riemannian metric. In fact, it does not even correspond to a general Finslerian metric, which 
in general appear in premetric electrodynamics [17j.[18j. 



5 Special axion fields 

5.1 Axion field on Minkowski spacetime 

In this section, we restrict for simplicity to the Minkowski spacetime. In the Cartesian coordinates, 
g v J = diag(l, — 1, —1, —1) with g = —1, so the dispersion relation ([36]) takes the form 

q 4 + (V\ m V>' m ) q 2 ~ (VV9™) 2 = . (38) 
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Note that now all scalar products are taken with respect to the constant Minkowski metric. We apply 
the (1 + 3) splitting and denote 

q=(w,k), fc = |k|, (39) 

and 

ip ti = (fi, m) , m = |m| . (40) 

Denote by a the angle between the vectors m and k. For a complex vector q 1 , the usual analytic 
extension for a is assumed. Due to the symmetry of (|38[) under the reflection ipi — > —ipi, we can deal 
locally with the case fj, > 0. 

In the (1 + 3) notation, the dispersion relation (|3"5)) takes the form 

(w 2 - k 2 ) 2 + (/i 2 - to 2 ) (w 2 - k 2 ) - (wfj, + mk cos a) 2 = . (41) 
It is useful to express the dispersion relation in term of the phase velocity v p = w/k 

* ( 2 + w ) ~ 2Vp[ ^ cosa + i 1 i + w sin2 a ) = ■ (42) 

We observe now that, in the left hand side of (|42|) . the dependence on the angle a can be removed only 
if the vector m can be taken equal to zero. Such coordinates can be chosen if and only if the 4-covector 
if; i is timelike. In this case, the 5*0(3) invariance is preserved. Alternatively, for an arbitrary null or 
spacelike covector ipj, the rotational symmetry is violated. 

Additionally, let the term linear in v p cannot be removed, i.e., the parameters m or fi cannot be 
chosen be equal to zero. Consequently, v p and (—v p ) cannot satisfy simultaneously the same dispersion 
relation. It means that the future and the past light cones are not identical. Thus the time inver- 
sion symmetry (T- invariance) is violated. Since the whole equation (|38)) is PT-invariant, the parity 
invariance (P-invariance) is also violated. 

The axion field and the covector of its derivatives are assumed to change smoothly in the whole 
spacctime. Thus the spacetime itself is separated to distinct regions with different norms of the covector 
ip m . In every specific region, special coordinates can be chosen in order to simplify the parametrization 
of the covector field tp m . 



5.2 Axion field with a timelike derivative 

This model and its physics consequences was studied in the original version of axion modified electro- 
dynamics - the Carroll-Field- Jackiw model [2]. Consider a spacetime region, where the derivatives of 
the axion field compose a timelike covector 

ip, t V > . (43) 

Choose in this region the time coordinate axis to be directed along the covector tj)^. Consequently, this 
covector is parameterized now as 

ti = (ji, 0,0,0). (44) 

Due to the symmetries of the dispersion relation (|36|) . we can, without lose of generality, require /i > 0. 
In the original CFJ-model, the axion field was given as ip = fit with a constant parameter [i. Observe, 
however, that only the first order derivatives of the axion field are involved in the dispersion relation 
(|36p . So, in fact, we can deal with a more general case, where /i is an arbitrary function of a point. 
Substituting (|HJ into (|JT]) we get 

(w 2 - k 2 - fik)(w 2 - k 2 + fik) = 0. (45) 

The solutions of this equation are 



Ww= Tjk 2 +(ik, (2) w = -^k 2 + iik, 

= yjk 2 -fik, (i) W = -yjk 2 - flk. 



(46) 
(47) 
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Thus, for k > /i, we have four distinct real solutions: two positive and two negative. For k — fi, there 
are two real solutions of opposite signs and one double solution equal to zero. For k < /i, two solutions 
are real and two are pure imaginary. 

Geometrically, the solutions (|46l |47|) define two distinct double hypersurfaces 



and 




Fig. 1 The graphs represent two hypersurfaces corresponding to CFJ modified electrodynamics 

We depict on Fig. 1 these hypersurfaces in the coordinates (fci , k%, w) (the third coordinate k% is 
suppressed, as usual). The first picture corresponds to the solution (|4"o| . It is topologically equivalent 
to the ordinary light cone. In particular, two cones are joined by a unique point (w = 0, k = 0). When 
this point is removed the interior region bounded by the hypersurface is separated to two disjoint 
parts - the future and the past spacetime regions. The second hypersurface corresponding to (|T7|) is 
topologically different from the ordinary light cone structure. It consists of two pieces which are joined 
by a sphere (fc = jj,, w = 0) and by a point (k = 0,w — 0). Consequently the interior region cannot 
be separated into two disjoint parts even if the origin is removed. It means that the future (the upper 
region) and the past (the downer region) always connected and cannot be disjoint. In other words, the 
causality on this branch is violated. 

Since the equations (|48|49|) involve a term linear in the covector k the corresponding light hyper- 
surfaces cannot be associated with some pseudo-Riemannian optical metrics. 

The phase velocities of the waves are given by the expressions 

- f^, "\ - f^, (50) 
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which coincide with the corresponding formulas of [2]. One of the phase velocities is greater than the 
speed of light in a vacuum. It increases monotonically when the factor ji/k increases, i.e., when the 
parameter k tends to zero. The second phase velocity is less than the speed of light in a vacuum and 
monotonically decreases to zero when the parameter k tends to zero. 

Well known that the phase velocity does not completely characterize the energy propagation. An- 
other useful characteristic is the group velocity which is usually thought of as the velocity at which 
energy is propagated along a wave. This quantity is defined by the derivative 



dw 
~dk 



(51) 



From (1501) we have 

^== , t2j %= •;,„'' \ ■ (52) 



mi k + n/2 (2) k-/J,/2 

\Jk 2 + [ik \Jk 2 — fik 



Consequently, both group velocities are superluminal and monotonically increase when the parameter 
k tends to zero. Observe that on the second branch, both velocities are defined only for < fi/k < 1. 
For k fi, the phase velocity tends to zero while the group velocity goes to infinity. This behavior 
corresponds to transmission of energy with infinite velocity, i.e., indicates the runaway modes, see [2] 
and P3J. 

Consequently, the waves with k > /i propagate along two distinct future light hypersurfaces. This 
behavior corresponds to the birefringence effect known from classical optics. The topological type of 
the light hypersurfaces is different, however, from the ordinary light cones. For k < /i, the runaway 
modes emerge. 



5.3 Axion field with a spacelike derivative 

Let us consider a spacetime region where the axion field has a spacelike covector of derivative 

i>,mt m < • (53) 
By transformation of the coordinates, we can choose in the whole region a parametrization 

^=(0,771,0,0). (54) 

Also here the parameter m can be considered as a function of a point. We can restrict to m > 0. The 
dispersion relation (|41~j) takes now the form 



w 2 = k 2 H ± \l h m 2 k 2 cos 2 a . (55) 

2 V 4 

Observe that an inequality w 2 > k 2 > holds for (|55|) . Hence this equation for w has four real solutions 
for every values of parameters. Consequently the runaway solutions are absent. 
The corresponding light hypersurfaces are given by 

^ 2 -^-^-^-yfl + ^+§j=0. (56) 

and 

w 2 -k 2 -k 2 -kl-^^-f7^=0. (57) 

These hypersurfaces are depicted on Fig. 2. The first structure is topologically different from the 
standard one. Indeed, ([56]) has not real solutions for w = 0. Consequently, two branch given by (1561) 
are represented by two disjoint surfaces. In other words there is not, for this lightlike hypersurfaces, 
a way from the past to the future. The second equation (|5"7| , for w = 0, has a unique solution 
k\ = ki = &3 = 0. Thus we have on this branch the standard light cone topology with one point 
joining the past and the future spacetime regions. However, an ordinary definition of the causality is 
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Fig. 2 The graphs represent two double light hypersurlaces of the spacelike axion modified electrodynamics. 



not applicable in this case. Indeed, due to the inhomogeneity of the dispesion relation, a type of a light 
trajectory depends on the parametrization. The expressions (|56[ I57p contain non-polynomial terms. 
Thus also in this case, the optical metrics cannot be represented in a pseudo-Riemannian form. 
In term of the phase velocity, (|42|) is rewritten as 



777 

J 2 



Thus the phase velocities are expressed as 



1 



k 2 



0. 



(58) 



\ 




and 



(2), 



\ 



1 



2k 2 



k 2 \Ak 2 



(59) 



(60) 



These expressions depend explicitly on the angle a, so the 50(3) invariance is violated. However, due 
to the fact that the expression is invariant under the change v p — > —v p , the future and the past cone 
are the same. Thus the T-invariance and consequently the P-invariance are preserved. For small k 
the first phase velocity (|59[) goes to infinity while the second phase velocity (|60[) tends to zero. In the 
transversal direction, a = 7r/2, one of the phase velocities is greater and one is equal to the speed of 
light in a vacuum. 

The group velocities are expressed as 
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Both expressions tend to zero for small values of k. In transversal direction we have 

1 



v. 



(62) 

p 



i.e., one of the group velocities is less and one is equal to the speed of light in vacuum. 
5.4 Axion field with a lightlike derivative 

Consider a spacetime region, where the derivatives of the axion field compose a null covector 

V,mV> ,m = 0. (63) 
In this case, the general covariant dispersion relation (|36[) takes the form 

<Z 4 -0/W") 2 = 0. (64) 

This expression is readily factored as 

(q 2 -t m q m )(q 2 + tmq m )=0- (65) 

The expression in the left hand side does not have a defined sign for an arbitrary non zero covector q m . 
Thus the non-birefringence condition |16j is explicitly violated. Consequently, the birefringence effect 
emerges for arbitrary varying null axion fields. 
We can choose a parametrization 

V>,i = (to, to, 0, 0) to>0. (66) 

On a flat Minkowski manifold, the dispersion relation takes now the form 

(w 2 -/c 2 ) 2 -m 2 (w + fccosa) 2 = 0, (67) 

or, 

w 2 = k 2 ± m(w + k cos a) , (68) 
Consequently the light hypersurfaces are expressed as 

w 2 -kl-kl-kl+miw + kx) =0, (69) 

and 

w 2 - k\ - k\ - k\ - m(w + fcx) = 0. (70) 
The linear terms indicate that these expressions cannot be represented by pseudo-Riemannian metrics. 



The equations (|68|) have four solutions 



m m m2 to <i\ 771 m2 

y >w = h \i hK^+TOfccosa y >w = \ h k 2 + mk cos a , (71) 

2 V 4 2 V 4 



and 



®iu = h \ hfc 2 — to/ccosq; — \ h k 2 — mk cos a . (72) 

2 V 4 2 V 4 y ' 

These expressions are real for every values of the parameters. Consequently, for the axion field with 
lightlike derivative, the runaway solutions are absent. 

The light hypersurfaces are depicted on Fig. 3. The first picture corresponds to (|69[) . while the 
second one is for (|7D|) . The future and the past surfaces are contacted at points 

777, 777 

w = ±—, k x =±-, fc 2 = 0, k 3 = 0. (73) 
This topological structure is also different from the standard light cone structure. 
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Fig. 3 The graphs represent two light hypersurlaces of the lightlike axion modified electrodynamics 



The electromagnetic waves propagate with two different phase velocities 



in 



(1) ^i + Vi^ + T cosa+1 < (74) 



and 



() ^ = -2fc + Vi^"T cosa + 1 - (75) 

When the dimensionless parameter m/k increases, i.e., for small k, the first velocity (f?4")) goes to 
infinity, while the second one ([75)) tends to zero. 
The group velocities are expressed as 

m = mcosa + 2k _ (7g) 

V m 2 + 4km cos a + 4fc 2 

and 

o\ mcosa — 2/c . . 

v g = = . (77) 

V m 2 + 4fcm cos a — Ak 2 

For small k they tend to the same value cos a < 1. 



6 Conclusions 

We have considered a general phenomenological model of axion modified electrodynamics. It is shown 
that the axion modified electrodynamics can be treated as a special case of premetric electrodynamics. 
In this formalism, the axion field is not involved as an additional Chern-Simon term. Alternatively, it 
emerges as an irreducible part of a general constitutive tensor. We have derived a covariant dispersion 
relation of axion electrodynamics. For a varying axion field, it yields a modification of light cone. The 
birefringence effect indicates violation of the Lorentz invariance for timelike, spacelike and null covector 
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of axion field derivatives. This effect, however, completely different from the ordinary birefringence 
appearing in classical optics and in the premetric electrodynamics. It can be explicitly seen from the 
fact that the topological structure of the light hypersurface is different from the ordinary light cone 
structure. In addition, the optical metrics are not homogeneous in the wave covector so they are even 
non-Finslerian. 

Acknowledgements My thanks to Roman Jackiw, Friedrich Hehl and Volker Perlick for useful discussion. 
My deep acknowledgments to Yuri Obukhov for valuable corrections. 



Appendix A. Tensors and pseudotensors 

Since the axion field possesses a special transformational behavior, it is useful to have a precise meaning of all 
quantities in volv ed in the formalism. Although the notion of weighted tensors and pseudotensors is a classical 
subject, see [20] and also [7] for a modern treatment, some confusions in the basic definitions can be found in 
literature. We will characterize the geometrical quantities accordingly to the transformation properties of their 
components relative to transformations of coordinates. 
Let be given a smooth transformation 

x'^x'' =fV) (78) 

with a transformation matrix 

and its inverse L,;<\ Denote the determinant of the transformation matrix by J. It is expressed by 

J = det(r' l ) = l ei , J>kn ,e^ kL r\ L j ' j L k ' k L l \ , (80) 

where and tijhi are the permutation symbols. They take the constant values 0, ±1 in all coordinate systems. 
The relation (|80)) is equivalent to 

/i'k'l'j = ^Mjj'.jj'. L k' k L V i ^ ^ 

and 

ti'j'k'V — J^ijklLii 1 Lji 1 hy Ly . (82) 

Recall, see |20] . the definitions of the extended tensorial objects. 

Ordinary tensor has a set of components which are transformed as 

T*'">... -v V'i ■ ■ ■ L/ ■ ■ • T"'j... . (83) 

Tensor density of weight k is a set of components which are transformed as 

1 

J* 

So the ordinary tensors are tensor densities of zero weight 



T i - j ,...^- 1 -L i i ...L/...T i - j .... (84) 



Pseudotensor has a set of components which are transformed as 

r*'-,-,... - (sgnJ^V-. V-.-T%... (85) 
Pseudotensor density of weight k is a set of components which are transformed as 

r .... //.-•./.. ''.•./' (86) 



The ordinary pseudotensors are of zero weight. 
Let us start with an action 



/ Lvol, (87) 

J M 



A = 

which is an ordinary real number. Its integrand, Lvol, has to be a scalar valued invariant volume element, i.e., 
a p seudoscalar density. In the formalism of differential forms, it is an odd (twisted) 4-form, see for instance 
[14] . Here L is a scalar valued function (Lagrangian) while vol is a special invariant volume element defined 
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from the geometric quantities. On a 4D pseudo-Riemannian manifold, the standard invariant volume element 
is defined by the determinant of the metric tensor g = det(gij) 

(R^m) yol = ^/^ d 4 x ^ d A x = dx dx l dx 2 dx 3 _ (gg) 

In a pure tensorial form, it is equivalently rewritten as 

' fllem 'uoZ = y] — g y^eijkidx' dx 3 dx k dx'^j . (89) 

Due to (|81|) . the permutation symbol e jkl is a tensor density of weight (+1). Its "inverse" eiju is a tensor 
density of weight (—1). Consequently d 4 x is a scalar density of weight (—1). 
From the ordinary transformation law for the metric tensor 

9i'j> = gijL^'L 1 j, , (90) 
one readily has the transformation law for the determinant and its square root 



Thus g is a scalar density of weight (+2) while y/—g is a pseudoscalar density of weight (+1). 

One builds from these two objects a quantity lijhi = yj—gtijki, which is a (non- weighted) pseudotensor. 
Observe that the integrand of the action Lvol has a proper transformation behavior, it is a pseudoscalar of 
zero weight. Since the expression Ly/—g is the subject of variation, this pseudoscalar density of weight (+1) is 
often referred to as the Lagrangian density. 

Let us return now to the premetric electrodynamics action 



M 



A= / [Fij-tft+AiT d 4 *, (92) 



and, equivalently, 



A — I ( \x iihl FiiFu + ) d*x . (93) 



M 



2' 



Since d 4 x is a scalar density of weight (—1), the scalar integrand (the expression in the parenthesis) has to 
be treated as a pseudoscalar density of weight (+1). It is constructed from the ordinary tensors Fij and Ai 
and the pseudotensor densities of weight (+1) — TC 3 and J 1 . On a pseudo-Euclidean manifold, the ordinary 
tensors FP 3 and J 1 are extracted from them by TC 3 = H %3 yj~^g and J 1 = ,P yf^g. Consequently, x ljk \ which 
is defined by TC 3 = X Fkl, is a pseudotensor density of weight (+1). 
In axion electrodynamics, the constitutive tensor is given by 

(Max) ijkl ( „ ih n i l „ il r ,i k 1 / n X j/ic'^' 



X ={9 9 -9 9 ) V-9 + ■ (94) 

Since yj—g s a pseudoscalar density of weight (+1), the first term is a pseudotensor density of weight (+1). 
In the second term, e 3 is a tensor density of weight (+1). Consequently, ip is a pseudoscalar. In physics 
literature, it is called axion. 

The matrix M t3 — X ljkl( lk<li is a pseudotensor density of weight (+1). Its determinant is a scalar density of 
weight (+2). Consequently the adjoint matrix Aij is a pseudotensor density of weight (+1) while the function 
A(g) is a pseudoscalar density of weight (+1). The photon propagator derived for premetric electrodynamics 
in [15] is a pseudotensor of weight (—1). 
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